Decay of metastable current states in one-dimensional resonant tunneling devices 
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Current switching in a double-barrier resonant tunneling structure is studied in the regime where 
the current-voltage characteristic exhibits intrinsic bistability, so that in a certain range of bias two 
different steady states of current are possible. Near the upper boundary Vth of the bistable region 
the upper current state is metastable, and because of the shot noise it eventually decays to the 
stable lower current state. We find the time of this switching process in strip-shaped devices, with 
the width small compared to the length. As the bias V is tuned away from the boundary value 
Vth of the bistable region, the mean switching time r increases exponentially. We show that in 
long strips lnr oc (Vth ~ V) 5 ^, whereas in short strips lnr oc (Vth ~ V) 3 ^ 2 . The one-dimensional 
geometry of the problem enables us to obtain analytically exact expressions for both the exponential 
and the prefactor of r. Furthermore, we show that, depending on the parameters of the system, the 
switching can be initiated either inside the strip, or at its ends. 

PACS numbers: 73.40.Gk, 73. 21. Ac, 73.50.Td 



I. INTRODUCTION 

In the last two decades with advances of miniaturiza- 
tion techniques various resonant tunneling structures be- 
came the subject of intensive research. It was experimen- 
tally observed that the current-voltage characteristics of 
resonant tunneling devices, such as double-barrier reso- 
nant tunneling structurea 1 ! 2 ! 3 ! 4 ! 5 (DBRTS) and superlat- 
ticeS ) 6 i 7 i 8 i 9 exhibit intrinsic bistabilities. Namely, it was 
shown that for each value of bias in the bistable region 
of the I-V curve the current can take two different val- 
ues. A bistable I-V curve of a double-barrier structure 
derived theoretically in Ref. is depicted schematically 
in Fig. ^ It is theoretically established! 1 ^ that near 
the boundary V t h of the bistable region (e.g., at point 
A) the upper current state is metastable. In recent ex- 
periments^ the switching from the metastable to stable 
current state was studied in superlattices. In particular, 
the mean switching time r was measured. 

The problems of decay of metastable states were 
studied theoretically in various fields, such as con- 
densed matter physics ^iiiii^iiSiiLiLiS quantum field 
theoryj 20 i 21 i 22 i 23 i 24 and chemical kineticsi 25 ! 26 In the con- 
text of resonant tunneling structures it was addressed in 
Refs. ITTIIT^ A typical double-barrier resonant tunnel- 
ing structure consists of three semiconducting layers of 
GaAs separated by two insulating layers of GaAlAs. In 
the narrow middle layer of GaAs the electron motion in 
the direction normal to the layers is quantized, so that a 
quantum well is formed. In the bistable region (Fig. ^) 
the two current states correspond to two different values 
of electron density n in the well. The density n exhibits 
shot noise fluctuations caused by random events of tun- 
neling of electrons in and out of the well. Thus it becomes 
possible for the devic e to sw itch from one current state 
to the other. In Refs. I 111 1 2l it was shown that there are 
two regimes of current switching. In the case of rela- 
tively small samples the electrons spread uniformly over 



the well due to diffusion. Then the switching occurs si- 
multaneously in the entire area of the sample, and it was 
found that the switching time r is exponentially large, 
with lnr oc (Vth ~ V) 3 ^ 2 . In larger samples the density in 
the well is not uniform, and the switching occurs through 
nucleation mechanism: it initiates in a small region of the 
quantum well of the characteristic size tq, which then 
spreads rapidly to the entire sample. The size of the crit- 
ical nucleus of the stable current state was found to be 
bias dependent, tq oc (Vth — V) -1 / 4 . It has been shown 
that in the large-sample regime lnr oc (Vth ~ V)- 

These results were obtained for samples whose two lat- 
eral dimensions are comparable to each other. In this pa- 
per we study the interesting case of the devices of strip 
geometry with w <ti L, where w is the width and L is 
the length of the strip. In these devices there is a special 
regime where the applied bias is such that w <C r$ <C L. 
Since the scale tq gives the characteristic size of the den- 
sity fluctuations, in this case the density does not vary 




FIG. 1: Current-voltage characteristic of the DBRTS. The 
bistable region is present in the range of bias between Vth and 
Vth- The inset shows a sketch of generic potential u(n) in the 
bistable region. The points A and C on the I-V curve cor- 
respond to the local and global minima of u(n), respectively. 
The point B on the unstable current branch corresponds to 
the maximum of it. 
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across the strip, but only along it. Confining the switch- 
ing process to one dimension alters its properties signifi- 
cantly. Similarly to the two-dimensional case, the switch- 
ing time grows exponentially when bias is tuned inside 
the bistable region. However, the exponent follows a dif- 
ferent dependence which is found in Sec. IIIII 

We show that in long strips the nucleation can oc- 
cur either inside or at the ends of the device. It turns 
out that both nucleation regimes can be observed, but 
the respective decay times are dramatically different. To 
compare them, one needs to calculate the switching times 
with the prefactors. Unlike the case of two-dimensional 
structures, in one dimension the prefactors can be found 
exactly in the limit of long and narrow strip, Sec. IIVI 

The decay of the metastable current state is governed 
by shot noise, and therefore switching is a stochastic pro- 
cess. We describe this process using the Fokker-Planck 
equation approach. In Sec. llV Al we introduce the Fokker- 
Planck equation for tunneling and in-plane diffusion in 
DBRTS. This equation is used to calculate the prefactor 
of the mean switching time in strip-shaped (Sees. IIV Bl 
and II V 01) and ring-shaped devices fSec. UV Dll . We sum- 
marize our results and discuss their experimental impli- 
cations in Sec. 



II. STATIONARY DISTRIBUTION FUNCTION 
OF ELECTRON DENSITY 



It was shown in Ref.[l2|that the stationary distribution 
function of electron density n(r) in the quantum well 
takes the form 

Po{n} = e- F{n \ F{n} = Jd 2 r [u{n) + ry(Vn) 2 ] . 

(1) 

Here the integration is over the cross-section of the well. 
In the case of uniform density F{n} is determined by the 
effective potential u(n), which describes the tunneling 
between the well and the leads. The gradient term in the 
functional F{n} accounts for the diffusion in the quan- 
tum well. The constant 7/ is positive and proportional 
to the in-plane conductivity, 12 r\ oc a. By suppressing 
gradients of n, this term favors the states with uniform 
density. 

The shape of the effective potential u(n) in the bistable 
region is schematically shown in the inset of Fig.^ It has 
the local and global minima at points A and C, respec- 
tively. In the case of uniform density n, these minima 
result in peaks of the distribution function Po{n}, which 
correspond to the upper and lower branches of the I-V 
curve. Point B on the unstable current branch corre- 
sponds to a maximum of u. 

In the vicinity of the threshold voltage Vtt the effective 
potential u(n) can be approximated by a cubic polyno- 
mial 

u(n) = -a(n-n t h) + -^(n-n t h) 3 , a = a(V t h-V). (2) 



Here 7 and a are positive constants. The voltage depen- 
dence of a ensures that the local minimum of it, which 
corresponds to the metastable state, disappears at the 
bistability threshold. The threshold density n t h is de- 
fined as the density at point V — Vth, where the lo- 
cal minimum A of u(n) disappears by merging with the 
maximum B. 

If the system is in the local minimum of u(n), it will 
eventually decay to the global minimum. In the limit 
of large conductivity a or small sample size the electron 
density n in the quantum well is uniform, and the gra- 
dient term in F{n}, Eq. (JIJ, can be omitted. In this 
case Eq. Q is simplified: the distribution function Pq 
is described by the only variable n and takes the form 
Po(n) = e~ Su ( n >, where S is the area of the sample. Ini- 
tially, the system is in the local minimum of u(n), see 
point A in the inset of Fig. In order to switch to the 
global minimum the system has to pass through point B. 
As it follows from the expression for Po(n), the probabil- 
ity of reaching point B is exponentially small, with the 
exponent determined by the barrier height S(ub — u-a)- 
The latter can be easily found from the expansion 
and the mean switching time tq takes the formic 



t = t exp 
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1/2 



(3) 



Here Tq is a preexponential factor, and we assumed that 
the cross-section of the sample has rectangular shape 
with the width w and the length L. 

Expression © is valid as long as L, w <C Tq, where 



1/4 



n> = 1 — oc {Vth - V) 
^7 



(4) 



has the meaning of the characteristic spatial scale of 
stochastic fluctuations of electron density.-— The scale r$ 
can be tuned by changing the bias V. If L, w 3> r the 
switching occurs according to the nucleation scenario. In 
this case the critical switching density is first achieved in 
a small part of the sample of size ~ r . After stochastic 
creation of the critical nucleus, it grows rapidly in size 
until it occupies the entire sample. 

In this paper we consider the case of a very narrow 
strip, w <C L. In the regime when the bias is such that 
w <C tq, the density may change only along the strip, 
and the problem becomes one-dimensional^ 

In the following it will be convenient to express the 
density n(x) in terms of a dimensionless function z(£) 
that vanishes at the minimum of u(n), namely, 



n(x) = n min - 2 x j—z{x/r ). 



(5) 



Here the density at the minimum n m j n = nth + \/ a /l- 
Substituting Eq. (J5J into the functional F in Eq. (QJ, we 
find 



(6) 
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where prime denotes differentiation with respect to the 
dimensionless coordinate £ = x/r along the strip. The 
characteristic value of the functional F is given by the 
parameter 



(a) 



Ui 



^3/4 



5/4 



(7) 



Its value depends on bias as U\ oc [Vth — V) 5 / 4 . 



III. THE EXPONENT OF THE MEAN 
SWITCHING TIME IN LONG STRIPS 

In this section we consider the regime ro ?S L. In this 
case, the density fluctuations along the strip result in a 
significant change of the mean switching time, and the 
result J3J is no longer applicable. The one-dimensional 
nature of this problem allows us to obtain an explicit 
expression for the exponent of the mean switching time 
r. 

We begin by finding the minimum and the saddle 
points of the functional F in Eq. © . They can be found 
from the condition SF/Sz = 0, that is, 



fz 



+ z 



0. 



(8) 



This equation should be solved with the boundary con- 
ditions z'(Q) — z'(A) = which account for the fact that 
there is no current flowing through the ends of the strip. 
Here we have introduced the dimensionless length of the 
strip A = L/r . 

Equation ijHJ can be interpreted as the equation of mo- 
tion of a classical particle with unit mass in the potential 
u(z) = -z 2 /2 + z 3 /3, see the inset of Fig. E^a). In this 
analogy the coordinate £ plays the role of time. There 
are two obvious solutions, z(£) = and z(£) = 1, cor- 
responding to the particle staying at the maximum and 
minimum of u(z), respectively. The minimum of F {z} is 
obviously given by z(£) — 0, since we defined z in such 
a way that z = at the minimum. The other solution, 
z(£) = 1, is a saddle point of F{z}. 

Apart from the two trivial solutions, equation |(HJ| may 
have ^-dependent solutions corresponding to a moving 
particle. The boundary conditions z'(0) — z'(A) = 
require zero velocity at the moments £ = and £ = A. 
Thus the particle performs oscillatory motion between 
turning points z(0) and z(A). If the particle starts at 
z(0) = cq in the range < cq < 3/2, one can easily find 
the other turning point c from the condition u(c) = u(cq). 
This equation has two solutions: 



3 - 2c ± >/-12cg-l-12co + 9 



The turning point z(A) corresponds to the positive root 
c+. 

The time A required for the particle to travel from one 
turning point to the other is obviously one half of the 
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FIG. 2: (a) Saddle point solution z 3 (£) for strips of length 
A = n, 3.2, 3.5, 4, 5, and oo. The allowed motion of a classical 
particle in the potential u(z) = — z 2 /2 + z 3 /3 is shown in 
the inset, (b) The dependence 1(A) defined by Eq. (TT31 . At 
A < 7T it is linear: 1(A) = A/6; while for A > n it is described 
by Eq. 11141 . At A — > oo the saddle point z s (£) is given by 
Eq. JTTJ, and I = 3/5. 



period of oscillations. The period of small amplitude os- 
cillations when the particle starts close to the minimum 
of the potential u(z), i.e. at cq — > 1, equals 2ir. The 
period monotonically grows to infinity as cq — * 3/2 or 0. 
Therefore, if A < tt there are no ^-dependent solutions, 
and z s (£) = 1. At A > 7T we have an additional sad- 
dle point z s (£) corresponding to a particle moving from 
z s (0) = c to z s (A) = c+. 

Equation (JSJ) can be solved analytically in terms of the 
elliptic integrals. In particular, the inverse function of 
z s (£) has the form 



V6 



F arcsin 



cq - z s 
Co - c+ ' 



c - c + 



Co 



(9) 

where F(ip, k) is the elliptic integral of the first kindi^ 
Using Eq. @ the length of the strip A can be expressed 
in terms of Cq as 



A = 



v/6 



-.F 



cq - c+ 
Co - c_ 



(10) 



Solving equation I|1U|) with respect to Co one can obtain 
the dependence co(A). Substituting it into Eq. and 
inverting £(z s ) one obtains the saddle point z s (£) for a 
given A. The saddle-point solutions z s (£) for several val- 
ues of A are shown in Fig. HJa) . 

For a strip of infinite length the boundary conditions 
take the form Co = z s (0) = 3/2 and c + = z s (oo) = 0. 
Then expression 10 can be significantly simplified and 



4 



yields 



2cosh 2 (£/2) 



(11) 



This solution can be easily verified by substitution into 
Eq. ©. 

The mean switching time r is given by 



r = r*e F ^\ 



(12) 



where r* is a preexponential factor. Using the expression 
for F{z}, Eq. iJBJ, we find 



F{z s } = I7i/(A), /(A) 



At A < 7r, the only saddle point is z s = 1, and therefore 
1(A) = A/6. At A > 7r the saddle point solution z s (Q 
given by Eq. @ corresponds to a smaller value of /. 
In particular, at A — > oo the solution i|ll|) gives / = 
3/5. In the intermediate region the integral 1(A) can be 
evaluated analytically, 




I 



A r 

60 ^ 
, \/6 



6(4 + c+ +c - c-) +9] 



,Jco-c-E 



cq - c + 
c - c_ 



(14) 



Here E(k) is the complete elliptic integral of the second 
kindiSS The dependence 1(A) is plotted in Fig.J^b). 

Note that apart from the saddle point described by 
Eq. © , for A > 2tt there are additional ^-dependent sad- 
dle points. For example, there is a ^-dependent solution 
corresponding to a particle moving from z s (0) = Co to 
z s (A) = c + and back to cq. This saddle point is respon- 
sible for the processes of switching inside the strip, see 
Sec. lIV (Jl At A > 37r another saddle point appears which 
corresponds to a particle moving from cq to c+ , returning 
to Co, and back to c + . In general, for A between irm and 
7r(m + 1) there exist m different ^-dependent solutions. 
However, the additional saddle points give larger values 
of F{z s } and thus do not affect the switching time. 



IV. PREFACTOR OF THE SWITCHING TIME 

The exponential dependence of the mean switching 
time (|12H was obtained from the stationary distribution 
function Po{n} of electron density. However, to calcu- 
late the prefactor r*, understanding of the time evolu- 
tion of the distribution function P{n(x),t\ is also re- 
quired. When electrons tunnel in or out of the well, the 
density n(x) changes in very small increments. Thus the 
dynamics of P{n(x),t} is described by a Fokker-Planck 
equation^ If the strip is very short, the density in the 
quantum well is uniform, and the system dynamics is 
described by a single variable n. In this case the Fokker- 
Planck equation for P(n, t) essentially coincides with the 



one for a small two-dimensional sample. Then by using 
conventional techniques^ the prefactor of the switching 
time can be found— as 



2vr 



(15) 



This result is correct as long as A <C 1 . 

In longer strips the density fluctuates, and therefore 
the Fokker-Planck equation with one variable n can- 
not adequately describe the evolution of the distribution 
function P. In this case P is a functional of n(x), and 
the Fokker-Planck equation is multidimensional. Then to 
study the decay of metastable current states one should 
use a more sophisticated method. 



Fokker-Planck approach to current switching in 
DBRTS 



To find the prefactor of the mean switching time for the 
system described by the multidimensional Fokker-Planck 
equation, we use the first passage time techniquei 30 i 31 
This method is applied to the Fokker-Planck equation in 
the form 



ap(x, t) 

di 



£P(x,t), 



(16a) 



£ = -E^*W+E*Sr^' xl - (1,i, '> 



h3 



dxidxj 



where the matrix Dij represents the generic diffusion co- 
efficient, and K is the drift field. Assuming that the 
system has a metastable state, one can consider its do- 
main of attraction Q with the domain boundary dfl being 
the separatrix of the field K. For the stochastic process 
described by Eq. Ijl6|l . the mean time of the first pas- 
sage out of the domain £1 has been found in Refs. I30l3ll 
The mean switching time is obtained as doubled mean 
first-passage time and takes the form^ 



r = — 



2 f n d d xP (x) 

E/ an ^EA,(x)P (x)^M 

i i 



(17) 



Here Po is the stationary solution of Eq. I|16f) . The form 
function /(x) is defined as a stationary solution of the 
adjoint equation, 



dxj 



= 0. 



£ f /(x,t) = E +E ^(x) 

(18) 

It satisfies the boundary conditions /(x) = 1 inside the 
domain O and f(x) = at the boundary dil. 

The Fokker-Planck equation for the case of non- 
uniform electron density n(r) in the well was obtained 
in Rcf. 12. It takes the form 

dP{n,t} 

M 



u'(n) 



2rjS7 2 n 



— )P{n,t}. 

on J 

(19) 
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It is easy to check that the stationary solution of this 
equation is given by Eq. (JIJ. 

Close to the threshold u(n) can be approximated by 
the expansion @, and the Fokker-Planck equation 119|) 
can be represented in terms of dimensionless variable 
z(£), Eq. JHJ), as 



dP{z,Q} 
56 



z - z > + ±-±y {z ,e } . 

(2°) 

where = 27rf/rQ is the dimensionless time. The sta- 
tionary solution Pq{z) of this equation is e~ F with F{z} 
given by Eq. ©■ Equation 1(201) has an infinite number 
of variables, since the density z(£) is different at every 
point. From now on we can consider a purely mathe- 
matical problem of decay of a metastable state for the 
system governed by dimensionless equation (|20|l . This 
equation is rather generic and we expect it to describe 
other one-dimensional problems of decay of metastable 
states. 



B. Nucleation at an end of a strip 

To find the prefactor of the mean switching time 1(12|) 
we use expression (|17|l . We evaluate both integrals in 
Eq. (|17J) in Gaussian approximation. The integral in the 
numerator of Eq. i|17|) is dominated by the minimum of 
F{z}. To find the expression for F in the vicinity of its 
minimum, it is convenient to use Fourier expansion 



(21) 

Substituting this expansion into Eq. ©, we find that up 
to quadratic in Xi terms, 



i=0 

where we defined A^ = (wi/A) 2 + 1. Using this expres- 
sion for F, one can easily evaluate the integral in the 
numerator of Eq. (|17|l , 



-_F(x) 



n 



2tt 



(23) 



The integral in the denominator of Eq. 1)17(1 is domi- 
nated by the saddle point z s (£) of the functional F{z}. 
To find the expression for F in the vicinity of the saddle 
point it is convenient to expand z(£) near z s (£) as 



K0 = *.(0 + 5><&(0. 



(24) 



i=0 



Here <pi (£) are the normalized solutions of the eigenvalue 
problem 



-_-2z s (£) + l)^(0 = Ai<M£), 



with the boundary conditions 



#(0)=0, #(A)=0. 



(26) 



Substituting Eq. J23J) into F{z} given by Eq. ©, and 
expanding near x = up to the second order in Xi, we 
find 



F{x} = £/i/(A) 



Ui 



-.2 



(27) 



j=0 



where the first term is given by Eq. 1(13(1 . 

Expression 1(27(1 implies that it is convenient to calcu- 
late the integral in the denominator of Eq. ((171) in terms 
of variables Xi. Since the expansion coefficients Xi are 
related to Xi by orthogonal trasformation, the integrals 
over Xi in the denominator of Eq. ((17|) can be replaced 
by those over ii. 

To calculate the integral in the denominator of Eq. I(17|) 
we also need to find and df/dxj. They can be ob- 
tained from the x-representation of the Fokker-Planck 
equation J3UJ. The adjoint to the operator C of Eq. lt2"U|) 
in the x-representation can be written as 



i=0 



OXi 



1 d 2 

Thdi 2 



(28) 



Here the terms of higher orders in Xi were neglected. 

In the denominator of Eq. ((17|) the boundary d£l is 
orthogonal to the unstable direction xq on the saddle. 
Therefore the sum over i reduces to the only term i = 0. 
Using the definition of D. L j we find from Eq. ((28|) that 
Dij — U^ l 8ij. The sum over j then reduces to a sin- 
gle term j = 0, and we need to find only the deriva- 
tive df/dio, which is given by equation ((18(1 in x- 
representation. The latter equation takes the form 



A x + — — \ 

00 Ui dxo J dxo 



0. 



(29) 



Since the saddle point is unstable in the io-direction, 
the eigenvalue Ao is negative. Thus equation l(29|) can be 
solved with the boundary conditions / = 1 at xq — > — oo 
and / = at the domain boundary x Q = required by 
the definition of /(x), cf. Eq. ((18(1 . As a result, we obtain 
df/dxo = ~{2\~X \Ui/tt) 1 / 2 at the saddle point. 

Substituting Eqs. (|23J) and |23 into Eq. (T^l we ex- 
press the switching time in a strip-shaped device r s as 



T _ T * e u im 




-T s , (30) 



where 




(31) 



Thus the evaluation of the prefactor of the switching time 
reduces to solving the eigenvalue problem 1(25(1 . ((26(1 . 
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1. Very long strip 



We first consider the limit of a very long strip, A ^> 1 . 
Then the saddle point z s (£) is given by Eq. so the 

equation l|25|) takes the form 



+ 1 mo = kMO, (32) 



d 2 3 
de ~ cosh 2 (£/2) 

The eigenvalue problem 1|32|) with boundary condi- 
tions 1(26(1 is solved analytically in the Appendix m In 
particular, we find that the discrete spectrum consists 
of two bound states with eigenvalues Ao = —5/4 and 
Ai = 3/4. At finite A the continuous spectrum is dis- 
cretized according to 



A; 



1 



1 



where the scattering phase shift S(q) has the form 

3 

6 (q) = arctan — . 

n— 1 



(33) 



(34) 



To evaluate T s it is convenient to take the logarithm 
of Eq. I|31[l . thereby converting the infinite product to a 
sum. In the limit A — > oo, we find 



InT, 



1 

7T ./II 



qS(q) 



1 



r 



dq 



ilnA. 
2 15 



(35) 



Using this result for T s and the fact that in a long strip 
Ai — > 1, we find the mean time of switching in a very 
long strip 



V2 . 
— r exp 



3Ui 
~5~ 



(36) 



It is important to note that the above calculation ac- 
counts for only one of two equivalent saddle points in a 
strip-shaped device. Indeed, the saddle point fluctuation 
|[n]| is positioned near the left end of the sample, £ = 0. 
Alternatively, the density fluctuation could occur at the 
right end of the device. This possibility is formally de- 
scribed by considering the saddle point z s (A — £). Since 
the two types of processes are equivalent, the total rate 
of switching at both ends of the device is 2/r e . 



2. Strip of arbitrary length 

Another regime in which analytical expression for r s 
can be obtained is that of A < it. In this case the saddle- 
point density is uniform, z s (£) = 1. Substituting it into 
equation (|2*5|l we find Ai = Ai — 2. Then expression 1(50) 1. 
(|31|l simplifies significantly, and the prefactor of t s be- 
comes 



n 



A 2 



(ni) 2 + A 2 ' 




FIG. 3: Prefactor r* of the mean switching time in units 
Tq vs. the dimensionless length of the strip A. At A < n the 
prefactor is given by Eq. 13711 . whereas above n it is calculated 
numerically. At A — > oo the prefactor approaches v2/5 in 
agreement with Eq. 13611 . 



Evaluating the infinite product we obtain 



sin A 
sinh A 



(37) 



This result generalizes the formula J2J) for the switching 
time in short devices, A <C 1, to the case of any A < tt. 

For a strip of finite length A > 7r, the saddle point den- 
sity z s (£) is given by Eqs. JOJ and fTPf. In this case the 
eigenvalues Ai are obtained by solving numerically the 
eigenvalue problem l|25|l with boundary conditions 
for a given A. Substituting them into expression 1(30(1 . 
(|31|l we find the switching time. The result for the pre- 
factor of Ts as a function of A is shown in Fig. 

Let us discuss the behavior of the prefactor at A = 
tt. At A < 7T the functional F{z} has only the uniform 
saddle point z s (£) = 1, whereas at A > tt there are two 
saddle points, the uniform one and the ^-dependent one, 
see Fig. EJa). These two saddle points merge at A = n. 
Formally, this gives rise to the fact that Ai vanishes at 
A = 7r, and therefore the prefactor r* = 0. We believe 
that the prefactor remains nonzero at A = n: however, 
the evaluation of the prefactor in this case requires a more 
careful treatment than the Gaussian approximation used 
in this paper. We leave this problem for future work. 

It is worth mentioning that the non-monotonic behav- 
ior of the prefactor, see Fig. does not result in non- 
monotonic dependence of the switching time on sample 
length A. The reason is that the prefactor r* is mul- 
tiplied by the very large exponential [see Eqs. ltT^)l - l(T5)) 
and Fig.|2fb)] which monotonically grows with A. There- 
fore, the expression (|12|l with the prefactor r* and expo- 
nent (|13|l is monotonic everywhere except for the very 
narrow region tt — 3/U± A < tt. As discussed in the 
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previous paragraph, a more accurate evaluation of the 
prefactor is required to evaluate the prefactor near A = n, 
and we expect that it would restore the monotonicity of 
t as the function of A. 



C. Nucleation inside long strips and in rings of 
large circumference 

Apart from the ends of a strip, nucleation can occur 
inside the sample. Such processes are most important in 
ring-shaped samples, which have no ends. Thus to study 
interior switching we model the sample by a strip with 
periodic boundary conditions. In the following it will be 
convenient to consider a strip of length 2 A with boundary 
conditions z(—A) = z(A), z'(—A) = z'(A). 

To obtain the switching time we first need to find the 
saddle point of F{z}. In the mechanical analogy used 
in Sec. IIIII [see also the inset of Fig. IHa)] the periodic 
boundary conditions for the saddle-point equation corre- 
spond to a full period of oscillations of the particle, rather 
than half period as in the case of nucleation at the ends 
in Sec. lIVBl One such solution is given by on the 

interval — A < £ < A, where z s (£) is defined by Eqs. @ 
and (JTUJ, see also Fig.^a). [Additional solutions are ob- 
tained by shifts z s (£) z s (£+A£).] To find the exponent 
of the switching time, we need to calculate F{z s } using 
Eq. jnj. Since the saddle-point solution 2 S (|£|) is a sym- 
metric function on the interval —A < £ < A, the integral 
in F{z s } is doubled compared to that of end switching. 
Thus the time of the switching inside the strip takes the 
form 



*2I(A)U 1 



r„ e 



(38) 



cf. Eq. fifty . 

The calculation of the prefactor r* is similar to the one 
for the switching at an end of the strip fSec. IIVB|l . As 
in Sec. IIV 51 it is convenient to evaluate the denominator 
of Eq. ||T2J by expanding z(£), Eq. where <fo(£) are 
again the eigenfunctions of equation but now with 
periodic boundary conditions 



M-A) = &(A), 4>'i(—A) = #(A). 



(39) 



The eigenvalue problem l|32|) . I|39|l is solved in the Ap- 
pendix D The discrete spectrum consists of one negative 
eigenvalue Ao = —5/4, one zero eigenvalue Ai = 0, and 
A2 = 3/4; all other eigenvalues are positive and belong to 
the quasicontinuous spectrum, Aj > 1 for i > 2. In the 
denominator of Eq. (|17|l the integration over the ampli- 
tudes X2,X3,X4, . . . of the modes with positive eigenval- 
ues is easily performed in Gaussian approximation. The 
integration over the amplitude x% of the zero mode is less 
trivial. 

Let us discuss the physical origin of the zero mode. The 
functional F{z} is translationally invariant on a ring, and 
therefore F{z s } does not change if the saddle-point solu- 
tion is shifted, z s (£) — > z s (£ + A£). In other words, the 



deformation 5z s = A£z' s (£) does not affect F. On the 
other hand, according to Eqs. I|24|) and (|27l) the deforma- 
tion Sz s = xi(j>i(£) of the saddle-point solution does not 
change F either, because Ai = 0. Since the zero mode is 
unique, we conclude 



A£4(O = zi0i(O- 



(40) 



Then it follows that^ z' s (£) = c 1 0i(£). The constant c x 
can be found from the normalization condition for the 
eigenfunctions 



(41) 



Using Eq. (|40|) and the relation z' s (£) — c 1 ^i(£), the in- 
tegral over the amplitude x± of the zero mode takes the 
form 



dx 1 = ci / d(A£) 



2c x A. 



(42) 



Expressions for df/dxo and Dij were obtained in 
Sec. IIVI independently of the exact form of the saddle- 
point density and are therefore still applicable. Substi- 
tuting expression l|23|) for the numerator of Eq. (|17|) and 
using Eq. (|42|l in the denominator, we find 



2|A |AiciAVc7r 




T 

-L r • 



(43) 



is similar to T s evaluated 
but with the eigenvalues A, and 



Here the infinite product T 
in Sec. ITVBl Eq. ffl} 

Ai calculated with periodic boundary conditions (|39|) . It 
is shown in the Appendix B that at A > 1 the quasicon- 
tinuous spectrum of the eigenvalue problem l|32l) , 1|39|) is 
still given by Eq. (|33|l . but becomes doubly degenerate. 
Therefore, in the limit of an infinite strip the product T r 
can be found as 



15' 



(44) 



cf. Eq. (|35|l . Also, substituting expression (|ll|l for z s (£) 
into Eq. (glj we find c x = ^6/5 at A > 1. Upon sub- 
stitution of these results into Eq. I)43|) the mean time of 
switching inside a long strip takes the form, 



Ti 



ibAy/u; 



exp 



6£/i 
"5" 



(45) 



where Ui and Tq are given respectively by Eqs. 10 and 
IjlSp. and A = L/2rQ. To obtain the exponent we used 
the fact that the integral 7 = 3/5 at A^ 1. 

In the case of a device of strip geometry, switching can 
be initiated both at the ends and inside the sample. The 
exponent of Ti is a factor of 2 larger than the exponent of 
T e [see Eq. (JSBJ], which makes the switching at the ends 
generally more favorable. On the other hand, the rate 
l/ri is proportional to the strip length L. Therefore in 
very long strips interior switching denominates. 
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D. Nucleation in rings of arbitrary circumference 



Expression l|45|l was obtained in the limit L > ro. In 
this section we discuss the case of a ring of finite circum- 
ference L < ro. We will consider separately the cases 
L < 2nr and L > 2nr . We will use the expression <|38|) 
for the exponential of the switching time, which was ob- 
tained for an arbitrary circumference of the ring. 

Let us start with the case L < Itttq. In dimensionless 
units it corresponds to the problem of a ring with the 
circumference 2A < 27r. Similarly to the case of a strip 
of length A < 7r, see Sec. lIVB"2l the saddle-point density 
is uniform, z s (£) — 1, and from equation l|25|l we obtain 
Xi = \i — 2. Analogously to the result (|30|) . 1311) for the 
strip geometry the prefactor takes the form 



(46) 



Since the boundary conditions for the ring are periodic, 
Eq. has two types of solutions, cos(7ri£/A) and 

sin(7ri£/A). The two solutions have the same eigenval- 
ues for any i > 0, and thus all egenvalues except for Ao 
are doubly degenerate. Repeating calculations similar to 
those for Eq. QJJ, we find 



0.04 




, sin A 
sinh A ' 



(47) 



This expression is similar to the result i|37p for the strip of 
length A < 7r. The absence of the square root in the right- 
hand side of Eq. (|4TJl is due to the double degeneracy of 
the eigenvalues. 

At A > 7r one needs to solve the eigenvalue prob- 
lem (|25ll with nonuniform saddle point £ s (£)j Eqs- @ 
and ljl Ojl , which can be done numerically. Then one sub- 
stitutes the eigenvalues A^ into Eq. I|43l) to find the pre- 
factor t*. As we found in the Appendix, in the limit 
A — * oo the three lowest eigenvalues Xi are nondegener- 
ate whereas the rest of them are doubly degenerate, be- 
cause the potential — 2z s (£) is then reflectionless. Unex- 
pectedly, our numerical calculation shows that the same 
property holds for any A > tt. Then similar to Eq. 144p 
we find T r = T^, where T s defined by Eq. I|3ip was com- 
puted in Sec. IIV B 21 The prefactor r* above A = tt is 
plotted as a function of A in Fig.^J One can see in Fig.|2| 
that at A — * oo the prefactor approaches as 1/A in 
agreement with Eq. I)45p . 

Contrary to the case of a strip, see Fig. [3 the prefactor 
r* does not vanish as A approaches tt from above. This is 
due to the fact that at A — * tt small A2 in the numerator of 
Eq. I|43l) is compensated by small c\ in the denominator. 
To see this let us consider the saddle point solution close 
toA = 7r. AtA = 7r + e with f « 1 this saddle point can 
be expanded up to its second harmonic as 

z s = 1 + k + ki cos(tt£/A) + k 2 cos(2tt£/A), (48) 




A 



FIG. 4: Prefactor r* of the mean switching time in units 
T o/VUi vs - the dimensionless circumference of the ring 2A 
above A — n. The dotted line shows the asymptote y^/lSA 
of t * at A — > 00, see Eq. I|45|l . 



where kq — —k\/2, k 2 = and e = 57tk^/12. Then 

from Eq. (|41|) we obtain ci = y/rrKi, and calculating A2 
up to second order in K\ we find A2 = 5k^/3. As a 

result the ratio y/x^/ci in Eq. I|43p remains finite, and 
the prefactor at A — > tt + takes the form 



V5 



2\/6sinh7r \/U\ 



(49) 



Similar to the case of a finite strip, the singular behav- 
ior of t* near A = tt is an artifact of Gaussian approx- 
imation. We expect that a more careful treatment will 
show that the prefactor remains nonzero and continuous 
at A = tt. 



V. DISCUSSION 

We have studied the mean time r of switching from 
the metastable to the stable state in one-dimensional 
double-barrier resonant tunneling structures. We calcu- 
lated both the exponentials and prefactors of r for the 
strip and ring geometries of the sample. In this section 
we discuss the behavior of the mean switching time in a 
strip-shaped sample depending on the distance from the 
threshold and the structural parameters of DBRTS. 

In the very vicinity of the threshold the exponent of the 
mean switching time is of order unity and the metastable 
state decays very rapidly. As the voltage V is tuned 
further inside the bistable region, the decay time becomes 
exponentially long. Our results are applicable in this 
regime. Thus the exponents in the expressions for the 
switching time in the regime of short strip, Eq. J2J, and in 
the regime of long strip, Eqs. (|36p and (|45p . must be much 
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(a) 



lnr A 



\ {v th -v?i* 




(V th - vfi 4. 



(b) 



V th V 



In r k 




FIG. 5: (a) Logarithm of the mean switching time r vs. volt- 
age at L < d. Near the threshold there is a region of 3/2- 
power law dependence of lnr (regime of short strip), then 
follows the region of 5/4-dependence corresponding to the 
switching at the ends in the regime of long strip, (b) Log- 
arithm of the switching time r vs. voltage at L d. Close 
to the threshold one first observes 5/4-power law behavior 
corresponding to the interior switching, then at voltage be- 
low V* follows the region of 5/4-dependence corresponding to 
switching at the ends. 



greater than unity. To check whether these conditions are 
satisfied, it is convenient to introduce a bias-independent 
characteristic length 



A- ( Wl 



1/5 



24 wrf 



1/5 



(50) 



chosen in such a way that the exponent in Eq. I|36(l takes 
the simple form (d/ro) 5 . 

The exponent of the mean switching time in the short- 
strip regime [Eq. (J3J] can be expressed in terms of d and 
tq as 



4Lwa 3 / 2 
37V2 



5 L / d 



(51) 



This exponent is much greater than unity at r <C 
(Ld 5 ) 1 / 6 . Since the regime of short strip is defined by 
ro L, it is present only if L <C d. In this case, 
starting at voltage difference Vth — V corresponding to 
ro ~ (id 5 ) 1 / 6 , there is a region of 3/2-power law de- 
pendence of lnr corresponding to the short-strip regime. 
Then, as Vth — V reaches the value corresponding to 
ro ~ L, follows the region of 5/4-dependence of lnr for 
the long-strip regime, see Fig.[5Ja). 

In the latter case the switching can be initiated either 
inside or at each of the two ends of the sample. Summing 



the rates of these processes we find 
1 2 1 



(52) 



To determine the type of 5/4-dependence that can be 
observed in the regime of long strip, ro <C L, one should 
compare the rates of switching inside and at the ends of 
the strip. Using Eqs. H36[) and l|45|) . the ratio of the rates 
r^ 1 and r" 1 can be expressed as 



exp 



ro J 2 r 



(53) 



Using this expression we conclude that at L <C d the 
switching always initiates at the ends rather than inside 
the strip. Thus the region of 5/4-dependence is described 
by Eq. ®. 

At L 3> d there is no region corresponding to the 
regime of a short strip. One can only observe 5/4-power 
law dependence of lnr corresponding to the long-strip 
regime. Analyzing Eq. I|53|l one can see that at L ^> d 
there are two distinct regions of 5/4-power law depen- 
dences. At 7'o d and very large L the switching initi- 
ates inside the strip. At very small ro the exponential in 
Eq. I|53|) becomes very large, and therefore the switching 
takes place at the ends. The voltage V* , at which the 
crossover between these two regions occurs, is given by 
the condition = r^ 1 applied to Eq. (|53() . 



V* w Vth 



v/./' V. h; d ) 



4/5 



(54) 



Thus we conclude that in the case of L ^> d, starting at 
V corresponding to r ~ d, one first observes 5/4-power 
law dependence of In r on voltage corresponding to the 
switching inside the strip, then at V below V* follows the 
region of 5/4-dependence corresponding to the switching 
at the ends, see Fig. Efb). 



Acknowledgments 

The authors are grateful to S. W. Teitsworth and M. B. 
Voloshin for fruitful discussions. O.A.T. acknowledges 
the hospitality of Argonne National Laboratory where 
part of this work was performed. This work was sup- 
ported by NSF Grant No. DMR-0214149 and the U.S. 
DOE, Office of Science, under Contract No. W-31-109- 
ENG-38. 



APPENDIX: EIGENFUNCTIONS AND 
EIGENVALUES OF EQUATION (l32l) 

In this Appendix we find the eigenfunctions of discrete 
spectrum and the phase shifts of continuous spectrum of 
equation (|32fl , as well as the respective eigenvalues A; . In 
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order to solve equation t'>2\) one needs to diagonalize the 
Hamiltonian 



H„ 



d 2 m(m + 1) 
W ~ 4cosh 2 (£/2) : 



with m = 3. 

The eigenstates of the Hamiltonian IjA.ljl can be ob- 
tained using algebraic techni que called supersymmetric 
quantum mechanics, see Refs. I33l34l35l We begin by in- 
troducing the raising and lowering operators 



a ± = T 



m £ 
— tanh — . 
2 2 



(A.2) 



Hamiltonians H m and H m -\ can be expressed in terms 
of a + and a~ as 



(m)^(m) 



Conversely, by acting with a~ on H m ^ 
with i > 0, one finds that a~\E'^ m ' ) are eigenf unctions of 
Hamiltonian if m _i with eigenvalues e| m . Therefore, we 

(A.l) 

conclude that in addition to its ground-state eigenvalue 
6q given by 1|A.4|) . the spectrum of i? ro consists of all 
the eigenvalues of if m _i. 

Using Eqs. (|A.4|) and (|A.5|I one can derive the eigen- 
functions and eigenvalues of H m from those of H m _\ and 
vice versa. To illustrate this, let us consider the Hamil- 
tonian Hq = — (d/d£) 2 . Its eigenfunctions are e ±lq ^, and 

the spectrum is given by eg = q 2 , where q is an arbitrary 
real wave number. Then using Eq. (|A.5(I we can find the 
eigenstates of Hamiltonian H\ . The continuous spectrum 
is obviously e g , and its normalized eigenfunctions are 



H m — a a — 



H m _i = a a + — . 



(A.3a) 
(A.3b) 



Let us consider the function = 1/ cosh m (£/2). It 
satisfies the condition a _ 4 r (^) = 0, and therefore is 
an eigenfunction of the Hamiltonian H m with eigenvalue 
e = — m 2 /4. Since has no zeros, it describes the 

ground state of H m . Thus we find the ground-state wave 
function and energy 



^/2ir{q 2 + 1/4) V ^ 



1 tanh | ) e l9 « 

2 2 , 



(A.6) 

The ground-state wave function and eigenvalue of Hi are 
given by Eq. i|A.4|) with m = 1. 



Applying the same technique two more times we find 
the bound states of H 3 , 



'r(m+i) 



1 



20Fr(m) cosh m (£/2) : 



£n = 



4 ' 
(A .4) 

Applying the raising operator a + to if m _i\I>^ m ^ = 
£ ( m 1 )vj / ( m 1 ) we find that a" 1 " 1 !'!" 1 ^ are eigenfunctions 



of Hamiltonian H m with eigenvalues e 
propriate normalization we obtain 



(m-l) 



Upon ap- 



1 



4^2 cosh 3 (£/2)' 
, ; . , _ VT5 sinh(£/2) 



*H0 



2\/2cosh 3 (£/2)' 
V3 



£ (3) _ _ 



5^ 



\/2cosh§ 4\/2cosh 3 | 



4' 



.(3) 



(A.7) 
(A.8) 



- 2 -(A.9) 



* i (m) (0 = 



1 



.(m-l) 



+ ,T,( m - 1 ) 



m 2 /4 



(-1 



,( m ) _ ,( m - x ) 
fc i — fc i-i 



(A.5) 



The eigenfunctions of continuous spectrum of H3 are 
given by 



v /2 7 r(g2 + 9 /4)( (? 2 + l)( (? 2 + 1/ 4 ) 



h - tanh 

d £ 2 2J \ d£ 



tanh - I I 1 — tanh - 

2 J V rf£ 2 2 



(A.10) 



r 



Noting that H3 coincides up to a constant with the 
Hamiltonian of equation (|32|l , we conclude that fa (£) = 

\E'P' > (£) . Therefore the discrete spectrum of Eq. I|32|l has 
three eigenvalues, 



A 



5 
4' 



Ai = 0, A 2 



The continuous spectrum of Eq. (|32J) is X q = 1 



The asymptotics of its eigenfunctions ijA.lOfl can be ex- 
pressed in terms of the scattering phase shifts 5{q) in the 
following way 



£ -> ±00. 



(A.12) 



(ATI) 

Here the absence of the reflected wave illustrates the fact 
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that the Hamiltonian H m describes scattering in reflec- 
tionlcss potential^ Comparing Eqs. (|A.10|I and l|A.12|) 
one can see that the phase shifts are given by Eq. (PIJ. 

For the evaluation of the prefactor of the switching 
time we need to study a finite system. In Sec. IIV CI we 
consider the system of length 2A with periodic boundary 
conditions l|39[l • At A 3> 1 the eigenfunctions of discrete 
spectrum of Eq. il-il'l) with the boundary conditions l|39|) 
are given by Eq. (|A.11|I . The quasicontinuous spectrum 
of a large hnite system can be obtained by applying pe- 
riodic boundary conditions (|39[l to the eigenfunctions 4> q 
in the asymptotic form l|A.12l) . As a result we find the 
wave number quantization 



where I is an integer. Because the phase shift <|34[1 is an 
odd function of q, it follows from Eq. (|A. 13|> that = 
—qi. Therefore, the eigenfunctions with wave numbers 
qi and g_/ have the same eigenvalue A = 1 + qf. Thus 



the quasicontinuous spectrum of Eq. (|32|l with periodic 
boundary conditions (|39[l is doubly degenerate. The two 
real solutions with the same X q can be represented as 
even and odd combinations of <f> q and 4>- q . 

Transcendental equation (|A. can be solved using 
iteration procedure. At A 3> 1 it is sufficient to perform 
only the first iteration. Substituting the result into Xi = 
1 + qf we obtain Eq. (j22J • 

In Sec. IIVBI we consider the system of length A with 
the boundary conditions l|26|) . All solutions of this prob- 
lem can be obtained from the solutions of the problem 
with the boundary conditions (|39[) . Indeed, any even so- 
lution of the eigenvalue problem (|32JI with the boundary 
conditions (|39|) satisfies the conditions (|26[) . Therefore, 
from Eqs. (|A.7fl - l|A.9l) and IjATlf) one concludes that the 
discrete spectrum consists of two eigenvalues —5/4 and 
3/4. Because only even combinations (j> q (^) + (f> q (~0 = 
fiqiO ~ <t>-q(0 satisfy the conditions J2HJ), the quasicon- 
tinuous spectrum is nondegenerate and given by Eq. 133|l . 
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